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Abstract. We show that monomial Cremona transformations of degree d in 



can have inverses whose degree d! is quite large (for d > 2, d! 
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occurs), and that the full list of possible degrees d! for fixed d and n does not 
always form an interval. An easy method for inverting the maps is presented. 
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1. Introduction: The problem of inverse degrees 
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A Cremona transformation is a birational morphism ip : P n ---» P n , where we fix 
n > 2 to avoid trivialities and sometimes write P£ to emphasize the field k. There 
are homogeneous polynomials gi of the same degree d in variables xq,x±, . . . ,x n 
such that if induces the partial function (xq : • • • : x n ) — > (go : • • • : g n ) ■ If 
nonconstant factors common to all the gi are canceled, the representation of ip is 
unique up to a nonzero constant factor, thus giving a well-defined degree d, even 
under independent coordinate changes in the domain and codomain. 

The inverse of ip also has a degree, here called the inverse degree d' of ip, so 
d" = d. The notation d! emphasizes that we are interested in which values of 
inverse degrees occur as <p ranges over Cremona transformations of degree d, in 
some fixed P£. In what follows it can be assumed that k (after extending) is 
algebraically closed. Classical tools, in a setting described for example in [12] . 
but going back to the earliest investigations by Cremona, Nother and Cayley, 
easily yield a bound d' < d n ~ l from the study of intersections with certain sets of 
hyperplanes. In P 2 this gives d! = d. In P 3 it gives y/d < d' < d 2 , a fact Cremona 
regarded as evident — see p. 278, 1. 6 of |3J. Cremona [3], [I] showed that all 
such d' occur here when d — 2 or d — 3. This is now established for all d in Pan 
[IT] , with examples that modify ones in [10] that were found to be flawed. In P n , 
^ . n > 3, many related questions remain open. 

We focus almost entirely on the much simpler monomial Cremona transforma- 
tions, where the defining polynomials g\ are monomial. These have been studied 
actively during the last decade — see for example Costa and Simis [2], Simis 
and Villarreal [H], and their references. As described in [S], with fuller details 
in Theorem 2.2 of [H], the inverse of any monomial Cremona transformation is 
monomial. Even in this special case, little is known about which inverse degrees 
are possible. Our results dash hopes for simple answers. 

A recent preprint [2] of Costa and Simis, based on [1], provides some definitive 
results on rational maps defined by quadratic monomials. One can see when such 
a map is Cremona, and if so quickly read off the degree of its inverse, by examining 
a graph whose vertices correspond to the variables, with the same number of edges 
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(allowing loops) which correspond to the quadratic monomials. Note that some 
articles such as [2] use n variables where we, following others, use n + 1. As a 
trivial application of the main result of [2] , one can see that a monomial Cremona 
transformation of degree 2 in P n , where n > 2, has inverse degree at most n, with 
the bound attained only when the graph is either a triangle, giving the classic 
Cremona transformation in P 2 , or a linear tree (chain) with a loop added at one of 
the two extremities, giving a map if : (xq : • • • : x n ) —> (xq : XqX\ : • • • : x n -\X n ). 
In particular, when n = 3 and d = 2, the maximum value of d' is 3 for monomial 
Cremona transformations, whereas without monomiality the maximum is 4. 

An example will be provided where the monomial Cremona transformations 
have inverse degrees fairly close to the upper bound of cf 1-1 , which in this special 
case ceases to be optimal for n > 2. Except for very small values of n and d, it 
is not practical to examine all monomial Cremona transformations to see which 
inverse degrees d! actually occur. Still, the few results observed from programs 
reveal that the d' do not always form a sequence: gaps can appear. This and 
other aspects relating to computation are discussed in the last section. 

2. Methods for obtaining inverses 

To prepare the ground for a formula for calculating inverse degrees, simplifying 
a method presented in Th. 2.2 of [14] , we establish notation and briefly explain 
easy ways to establish some of the most relevant results. 

The notation used to define ip will be changed slightly to admit m+1 monomials 
go, ■ ■ ■ ,g m m n+1 variables, so <p becomes a rational map from P n to the closure of 
its image in P m . In many sources monomials are taken to be monic, but we shall 
temporarily allow constant factors, to indicate how little effect this generalization 
has on the maps. Only the main case m = n will be described, roughly following 
[8]. The group of monomial Cremona transformations in P£, k an arbitrary field, 
is a split extension. First, the maps (x : • • • : x n ) H- (c\Xq : • • • : c n x n ), Ci G k*, 
where constant factors act trivially, form a normal subgroup T, a torus (k*) n . 
Multiplying by elements of T does not affect degrees of maps. In the group of 
monomial Cremona transformations, the monic monomials, which are the real 
focus of interest, form a complementary subgroup to T. As determined in [5], 
or from discussions about lattices below, this subgroup can be identified with 
Aut(Z n ), or GL n (Z). It acts on T as the full group of rational automorphisms. 

Henceforth, monomials will be monic. Given systems of coordinates, every 
rational map (f : P n --- » P m defined by monomials Q{ can be represented by its 
fi + lxm + 1 log-matrix A, whose j th column lists the exponents of the monomial 
Qj. A log-matrix A is stochastic or, for emphasis, ci-stochastic: every column has 
the same sum, which is d. Since only exponents of the gj are used, the field k is 
now irrelevant. The rational map p is unchanged if, in its log-matrix, all entries 
in a row are adjusted by the same amount, remaining in N. Thus it can usually 
be assumed that the minimal entry in each row of A is 0, in which case its column 
sums give the degree of <p, and we say that A is reduced. 

For greater flexibility, matrix entries will now be allowed to range over Z, and 
the <7i may be Laurent monomials. For each index i, let Ri be the 1-stochastic 
matrix whose entries in the i th row are 1 and whose other entries are 0. Stochastic 
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matrices will be called equivalent if they have the same size and their difference 
is an integral linear combination of the R4 or, in other words, if they induce 
the same rational map. Matrices are supposed be n + 1 x m + 1, except when 
monomial rational maps are composed, which corresponds to multiplying matrices 
of compatible sizes, then usually passing to the reduced matrix. 

To say A is <i-stochastic means 1A = dl, where 1 denotes an all 1 row vector, 
here of size n + 1. Then multiplying by A on the left induces a map of Z- 
lattices (free abelian groups) A™ — > Aq, where A™ denotes the sublattice of Z m+1 
consisting of column vectors whose entries sum to 0. Stochastic matrices are 
equivalent precisely when they induce the same lattice map. 

There are many criteria for deciding if a map ip is birational. Some that are 
specific to the monomial situation appear in [T3] and [T5]. If the above lattice 
map is not an isomorphism, the associated monomial map tp has no inverse, not 
even among rational maps. A proof can be extracted from Sections 1 and 2 of 
[13], but uses more machinery than would be expected for such a basic fact, so 
the following approach seems worth recording. The key result is: 

Theorem 1. If a birational map ip : P n — -> W , W closed in ¥ m , is defined by 
monomials, the induced lattice map A™ — > Aq must be surjective. 

Proof. By assumption, ip has a rational inverse ip, which must define an injective 
function on some nonempty open subset V of W via n + 1 homogeneous poly- 
nomials in variables yo, . . . y m coordinatizing P m , with open image U C P n on 
which ip is defined. All sets here are irreducible. If desired, one can assume that 
the coordinate functions Xi or yj vanish nowhere on U or V, to work with rings 
Ru, Rv of functions on U, V, generated by Laurent monomials of degree in the 
coordinates. Recall that Aq contains the image of the map on A™ induced by the 
log-matrix of ip. A suitable change of basis for Aq, and a shift to multiplicative 
language, produce Laurent monomials Zi that are free generators of Ru, such that 
certain powers zf* generate the same multiplicative group as that generated by 
the monomials ip*{yj). Although ip induces an isomorphism between k(U) and 
k(V) (quotient fields of Ru and Rv), the images of the already generate k(V). 
This forces \di\ = 1 for all i, so the above lattice map is surjective. □ 

A rational map ip : P n — » P n defined by monomials is birational precisely when 
it induces an automorphism of the lattice Aq, and the inverse must be monomial. 
Such maps have been dealt with before, most notably in [H]. Our main aim is 
to provide practical methods for obtaining inverse maps and their degrees, where 
the matrices can have entries in Z. 

Let ip be birational, with <i-stochastic log- matrix A acting on Z n+1 . The re- 
striction to Aq has determinant ±1 and, for every column vector v, Aw — dv lies 
in Aq. Thus det(A) = ±0? (cf. Lemma 1.2 of [13]). Now assume d ^ 0. Since A^ 1 
acts on Aq, the difference of any two columns of A~ l is integral (cf. Lemma 2.1 of 
[H]). This inverse is of the form d~ 1 A*, where A* has entries in Z. In each row 
i of A~ l , let 7*i be the least element, and define ki = —dri. All k{ lie in N, since 
A~ 1 A = I and n + 1 > 2. Using again the matrices Ri with i th row 1 to adjust 
rows, define B = (I+J^i k^R^A" 1 . Then B is the reduced matrix that represents 
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the inverse of f, so its column sums give the inverse degree. This yields a formula 
that may have theoretical use: 

Theorem 2. Let A be a d- stochastic matrix obtained from Laurent monomials 
that define a Cremona transformation. Then the inverse map is also monomial 
and, assuming d ^ 0, its degree is d' 1 — J2 r i> where denotes the least value 
occurring in the i th row of A -1 . 

To clarify this result, and to eliminate the apparent need to work with fractions, 
one should work within the framework of Laurent monomials in affine coordinates 
Xi = Xi/xQ. We present this indirectly via an easily applicable matrix-based 
approach, in which d and d' are regarded as functions on the group G of auto- 
morphisms of Aq, with d'(g) = d(g~ l ) G Z + . 

Using the ordered basis e$ — eo (0 < % < n + 1), G is identified with the matrix 
group GL n (Z). Starting with the matrix of any g G G, one prepends a row 
and column in the zero-th position, so that all entries of this column, and all 
column sums, are zero. Adjusting rows as before yields a reduced matrix A = A g 
which, as should be clear from earlier discussions, is the log-matrix of a monomial 
Cremona transformation if whose column sums give the degree d(f). Thus d(g) 
is the sum of n+ 1 nonnegative terms, one for each row of A, using the negative of 
the least term in the row, or if no negative terms are present. Similar formulas 
appear in Remark 3.2 of [8] in a loosely related context. 

Example 1. We pass back and forth between some matrices in GL^Z) and the 
related log- matrices for monomial Cremona transformations of P 2 (so d' — d). 
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Permutations of coordinate vectors induce a subgroup S'n+i of G = GL„(Z), 
and d is constant on each double coset S n+ igS n+ i, since this involves permuting 
rows and columns of stochastic matrices. By considering A gi A g2 , one sees that 
digig2) < digi)dig2). Also, for g S n+ ±, dig) > 2, and d(—I) = n. Such ideas, 
developed further, may well yield useful information. 



INVERSES OF MONOMIAL CREMONA TRANSFORMATIONS 



5 



3. Examples with large inverse degrees 

There are in the literature many examples of monomial Cremona transforma- 
tions with diverse features of interest — see for example Section 5 of [15J. An 
attempt to construct maps in P n of degree d with relatively large inverse degrees 
d' led to the following examples, which may be those with the largest possible d' . 
Extensive tests produced no other example (up to permutations) that attained 
or exceeded the value given in this section. We will proceed as in Theorem 2. 
This should be compared with the later approach. 

Example 2. For d,n>2, define <p : (xo : ■ ■ • : x n ) — > (xq : XqX\ : ■ ■ ■ : x^z\x n ), 
which has an upper-triangular log-matrix A whose nonzero entries are an = d, 
and, for i > 1, aj_i 5 j = d — 1, an — 1. Let B be the upper-triangular matrix 
with entries by = (1 — d)- 7 " 4 for i < j. It is easy to verify that dividing the first 
line of B by d produces A' 1 . For convenience, write c = d — 1. To illustrate, 
the matrices A and B (almost inverses) are shown below when n is 2 or 3. The 
degree of (p is d and the above criterion will, after simplification, give the inverse 
degree d' = 1 + c H h c n_1 , which is ^ ^ d > 2. 

For n of the form 2m, the raw sum for d! is + ^rj~ +c 2m ~ 1 +2 YlT=o c 2l+1 +0, 
whose fractional parts yield an alternating sum that meshes with the last sum. 
The case n = 2m + 1, d! = + ™ x + 2 Y1T=q ° 2l+1 + 0, ^ s a l mos t identical. 





' c+1 


c 


" 




" 1 


— c 


c 2 " 


For n = 2, A = 





1 


c 


,5 = 





1 


— c 










1 










1 



For n = 3, A 



c+1 


c 











1 


c 











1 


c 











1 



5 



1 


— c 


c 2 


-c 3 





1 


—c 


c 2 








1 


— c 











1 



4. Computational aspects 

All monomial Cremona transformations of degree 5 in P 3 can be generated 
from the 367,290 combinations (unordered selections) of 4 out of 56 monomials 
in xq, . . . ,Xs, discarding those that fail to define a birational map of degree 5. 
Only 11,496 maps survive. By using combinations, the effect of permuting the 
columns (four coordinates in the image) has been factored out, so that only one 
out of each class of 24 maps is counted. The S4 that permutes rows acts on these 
classes, in orbits (often regular) consisting of maps with the same inverse degree. 
The frequency of occurrence of each d', but not the orbits, of these 11,496 maps, 
is shown below. The d' for the 48,042 monomial Cremona transformations of 
degree 3 in P 4 are also shown. In that case, four kinds of maps are square-free, 
with descriptions given in Prop. 5.5 of |15j . 
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Table 1. All monomial Cremona transformations with d = 5 in P 3 . 



Value of d'\ 
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Table 2. All monomial Cremona transformations with d = 3 in P 4 . 



Value of d'\ 


2 


3 


4 


5 
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7 
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Freq. of d'\ 


432 


8670 


14640 


10920 


5820 


3720 


1200 
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15 
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840 


360 





240 





120 



One sees from orbit sizes that, in both cases, the only examples giving the max- 
imum value of d' are those described in the previous section, up to permutations 
of rows and columns. More importantly, even for these small values of d and n, 
the list for d! contains gaps instead of forming an interval. Similar exhaustive 
calculations show that for monomial Cremona transformations in P 3 , the next 
case d = 6 has values of d' that do form an interval [3, 31], and here each d' from 
28 to 31 arises from an essentially unique map. When d — 7, d' can assume any 
value in [3, 43] except 39 or 40, while for d = 8, d! G [4, 57], excluding only 54. A 
double gap first appears for d = 10, where d' cannot be 84, 87 or 88. 

Moving to P 4 , it is unlikely that theory will help explain the results. For 
example, when d = 4, the d! lie in [2,40], with gaps at 32, 34-36, and 38-39. 
When d = 5, the d! lie in [3, 85], with gaps at 63, 70, 72, 74-75, 77-80, and 82-84. 

The calculations were performed in SINGULAR [6], with about 100 lines of 
code, available on request. On an ordinary computer, this took a few minutes, a 
few hours, or (for d = 5 in P 4 ) about two days. Results from a similar program 
that randomly generated large numbers of Cremona transformations with fixed 
d and n had already hinted strongly where the smaller gaps lie. 

Since this seems to be the first record of systematic calculations made in this 
area, it seems appropriate to mention what data can and cannot be expected to 
be obtainable. Tables for values of d up to 12 in P 3 have been produced, with 
13 expected to be the last, as the naive algorithm used becomes too slow. An 
implementation to count the number of monomial Cremona transformations in 
P 4 with given multidegrees [10| d\ = d, o^, d 3 = d', and giving data for a few 
values of d, should soon follow. 

The diversity and richness of examples obtained already provide a counter- 
part and challenge to theory. To go further, it seems worthwhile to use random 
sampling to assemble extensive collections of certain kinds of Cremona transfor- 
mations, not just monomial ones, then search for new phenomena. This needs 
to be done with care. We point out where some simple approaches run into 
difficulties. 

The idea used above for random generation of monomial Cremona transforma- 
tions quickly becomes useless as d and n grow, as d-stochastic matrices giving 
Cremona transformations are rare, and those giving large values of d' are even 
rarer. Also, the model used to generate columns (d approximately independent 
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and uniform selections of a position) is strongly biased against outcomes where 
the nonzero values are concentrated in few entries, the ones most likely to give a 
d! that is unusually large or small. By repeated sampling until manjo monomial 
Cremona transformations were obtained, the maximum d! observed started to 
fall well below the bound given by the examples of the previous section, blocking 
attempts to clarify where further gaps might lie. To show just one small exam- 
ple, by generating only 10,000 Cremona transformations with d = 6 in P 3 we did 
come fairly close to the maximum value d' = 31 but, as typically happened, failed 
to observe the minimum, here d! = 3, although values just above this appeared 
frequently. 

Table 3. 10,000 samples (from 3,226,875 attempts) with d = 6 in P 3 . 



Value of d'\ 


4 
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10 


11 


Freq. of d': 


260 


270 


797 


838 


1169 


1737 


1031 


1342 
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19 
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224 
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158 
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36 


56 
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27 




64 


17 
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A better approach is to work in a fixed P n , using matrices in GL n (Z) to gen- 
erate monomial Cremona transformations with variable parameters (d, d') (or 
multidegrees) whose range of values should become clearer over time. A stream 
of matrices can be obtained from the identity matrix by making random choices 
to replace some line with that line plus a multiple of another line. Any matrix in 
GL n (Z) is reachable by this process: just consider how to row- reduce nonsingular 
matrices when working over the integers. For a slightly different explanation, see 
[TB] . In practice, best results were obtained using multiples restricted to have 
small absolute value (even 5 seems too large), which helped postpone arriving at 
a case with d or d' beyond the range considered, at which point the generation 
process was started afresh. This idea has the great advantage that it produces a 
usable example at each step. However, building a table of pairs (d, d') that occur 
in a certain range is a slow process, as the cases of most interest tend to appear 
rarely — in some cases, only once in more than 330 million trials. Several days 
of computation were enough to fill most of the possible positions repeatedly, but 
in no case was there compelling evidence to suggest whether some new apparent 
gap was real or illusory. 

Acknowledgment. A conversation with Ivan Pan helped clarify details 
about bounds for d' . He also sent an early version of [IT] . 
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